o 
o 



c^ 



density functionals 



o : Non-universality of commonly used correlation-energy 

u 
< 

CM. , , 

Jacob Katriel* Sudip Roy' and Michael Springborg+ 

^~-*' Physikalische und Theoretische Chemie, Universitat des Saarlandes 

S: 

D . 66123 Saarbriicken, Germany 

^ : February 2, 2008 

o 

^, 

■ Abstract 



The correlation energies of the hehum isoelectronic sequence and of Hooke's atom 
isoelectronic sequence have been evaluated using an assortment of local, gradient and 
meta- gradient density functionals. The results are compared with the exact correlation 
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vQ ■ energies, showing that while several of the more recent density functionals reproduce 



the exact correlation energies of the helium isoelectronic sequence rather closely, none 
is satisfactory for Hooke's atom isoelectronic sequence. It is argued that the uniformly 
^' acceptable results for the helium sequence can be explained through simple scaling 

Q-i' arguments that do not hold for Hooke's atom sequence, so that the latter system 

- ^_, , provides a more sensitive testing ground for approximate density functionals. This 

^ • state of affairs calls for further effort towards formulating correlation-energy density 



functionals that would be truly universal at least for spherically-symmetric two-fermion 
systems. 
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1 Introduction 

Density functional theory (dft) can be implemented in a variety of ways, all of which have 
their roots in the pre-Hohenberg-Kohn era. A fully consistent implementation requires the 
formulation of density functionals for the kinetic energy as well as for the exchange and 
correlation energies, none of which can be expressed in terms of the one-particle density 
in a straightforward way. The Kohn-Sham formulation depends on the construction of a 
local one-body potential that incorporates the effect of interelectronic repulsion, giving rise 
to a ground state density that coincides with that derived from the solution of the exact 
many-body hamiltonian. Finally, an approximate approach, based on the fact that the 
Hartree-Fock density is fairly close to the exact one, proposes that the exact ground state 
energy be evaluated by adding to the Hartree-Fock energy a correlation energy, calculated 
in terms of the Hartree-Fock density by means of an appropriate density functional. 

Since the ultimate correlation-energy density functional, whose existence is guaranteed 
by the Hohenberg-Kohn theorem, remains elusive, the development and testing of ever im- 
proved density functionals keeps being vigorously pursued. For an extensive recent study see 
Kurth, Perdew and Blaha fT| . In this vein, Chakravorty et al. '2] and Jarzecki and Davidson 
[3] compared the exact and dft correlation energies across the helium (and other) atomic iso- 
electronic sequences. In the former article the exact correlation energy for the ground state 
of the helium isoelectronic sequence (as well as for all other atomic isoelectronic sequences 
with up to 18 electrons) is compared with the values obtained using six density functionals 
(i.e., LYP i, CS 0, VWN-SPP |6; , PZ |7: , B88 , LC — the acronyms will be explained 
below). Of these, only the LYP correlation energies are uniformly close to the exact values 
over the range of nuclear charges examined. B88 grows to become twice as large as the exact 
correlation energy, and all the others seem to diverge for large Z. The latter article extends 
this study, considering P86 |10j, LYP \^^ (again), and PW91 [TT as well as various combi- 
nations of these with exchange functionals. While, as noted above, the correlation energy 
obtained with the LYP functional (evaluated for the Hartree-Fock density) follows the exact 
correlation energy rather closely, this is not the case if the same functional is combined with 



Becke's exchange functional. On the other hand, the combination BP86 (Becke's exchange 
with the P86 correlation functional) exhibits a rather close agreement with the exact corre- 
lation energies for the helium isoelectronic sequence. The overall conclusion of these authors 
is "all DFT functionals, [...] reveal problems in descriptions of the energy as a function of 
nuclear charge" [S]. Correlation energies were also evaluated by Umrigar and Gonze ^2] for 
several members of the helium isoelectronic sequence, in terms of both local and GGA-type 
density functionals. 

The high Z limit of the helium (and other) isoelectronic sequences within dft was carefully 
examined by Whittingham and Burke ^S], who pointed out that some correlation-energy 
functionals do not scale correctly to the high-density limit. In particular, the local-density 
approximation diverges logarithmically at large Z, and the Perdew-Wang (PW91) functional 
[TT] fails as well. They note that PBE [[Hj was designed to correct for this failure. 

The over-estimate of atomic correlation energies by local density approximations based 
on the homogeneous electron gas has been analyzed by Tong [15 , who attributed it to 
the continuous excitation spectrum of the latter, infinitely extending, system, that allows 
for a "soft" response to the internuclear repulsion. A modification of the local density 
approximation, in which an energy gap is introduced, was recently studied by Rey and Savin 
|16j . Their results suggest that the inhomogeneity needs to be introduced more explicitly. 

An analysis of the properties that correlation-energy density functionals should satisfy in 
order to present a valid behavior along isoelectronic sequences was presented by Staroverov et 
al. ^Tj, who, as Whittingham and Burke JH], studied the high Z limit of some two-electron 
atomic ions. Furthermore, an interesting methodological discussion of the construction of 
density functionals was recently presented by Perdew et al. ^H]- These authors make a 
strong case for nonempirical functionals and specifically recommend the PBE [l4,j and the 
TPSS [ini functionals. 

Hooke's atom consists of two "electrons" that repel one another like ordinary electrons 
but are bound to the "nucleus" via harmonic potentials. Its separability, in center-of-mass 
and relative coordinates, was pointed out by Kestner and Sinanoglu [20]! and the fact that 



for certain values of the force constant it is analytically solvable was noted by Kais et al. 
PT] and, more systematically, by Taut 22j. Several studies for Hooke's atom have been per- 
formed, aimed at deriving exact density functionals or Kohn-Sham type exchange-correlation 
potentials [23-25], most notably by Burke et al. [26-30] and Ludena et al. [31J. However, a 
systematic study over a broad range of force constants has not been presented, to the best 
of our knowledge, with the exception of the local-density calculations that were carried out 
by Huang and Umrigar [32] ■ 

In the present paper we examine an assortment of correlation-energy density functionals 
proposed over the last seventy(!) years. Virtually each one of them had been tested for the 
helium atom, and most of them had been tested for a range of atomic ions isoelectronic with 
helium. Needless to say, many of them have been applied, with varying degrees of success, 
to the study of heavier atoms as well as molecules, surfaces and solids, and share the credit 
for the immense impact of dft in all these areas. 

The Hohenberg-Kohn theorem implies that the ground-state correlation energy is a uni- 
versal functional of the ground state density. The examination of the extent to which the 
commonly used correlation-energy density functionals satisfy this property is the principal 
aim of the present paper. This is carried out by considering the performance of the presently 
selected set of correlation energy density functionals for Hooke's atom isoelectronic sequence. 
Since these density functionals were originally developed for Coulombic (atomic and molec- 
ular) systems, their performance for Hooke's atom provides a severe test of their universal 
validity. Moreover, Hooke's atom remains bound over the whole range, < fc < cxd, of the 
force constant. This fact allows correlation to be examined not only in the high k (weak 
correlation) limit but also in the low k (strong correlation) limit. The latter limit is not 
available for atomic systems, that become unbound below a certain critical nuclear charge. 
This property of Hooke's atom isoelectronic sequence is elaborated upon below. 

Another model system with similar properties is that of a two interacting particles con- 
fined to a spherical volume. Here, the radius of the sphere provides a parameter that, upon 
variation, can be used to change the system from the low correlation to the strong correla- 



tion limit. This system was recently studied in detail by Jung et al. [3^ who also compared 
the exact exchange and correlation energies with those obtained by some few approximate 
functionals. Thus, the combined results of their and our study provide detailed insight into 
the accuracies of currently applied approximate density functionals and, consequently, may 
be used in improving those. 

Ultimately, the physical reason behind the existence of correlation effects is the electron- 
electron interactions. Thus, the total correlation energy can be written in terms of a double- 
integral over two position coordinates, involving the exact and the Hartree-Fock first and 
second order reduced density matrices, 

Bo = //p(r-Op(r-.)f.(n.f.)dMd¥.. (1) 

where ec(n, ^2) is a unique correlation-energy density in a six- dimensional (ri, r2) space. On 
the other hand, in density functional theory it is common practice to write the correlation 
energy as 

E, = Jp{r)e,{r)d'f. (2) 

The transformation from Eq. (^ to Eq. Q is not unique. Therefore, different approximate 
expressions for edr) may look very different but nevertheless give very close total correlation 
energies E^ jM]- Thus, since the approximate correlation-energy functionals most often are 
given in terms of approximate forms for ^d'f^), it is not possible to compare those directly 
and thereby obtain information on their performance. 

Instead, only the total correlation energy is of physical relevance. Therefore, by consid- 
ering the total correlation energy for continuous classes of systems detailed information on 
the performance of approximate density functionals can be obtained, that ultimately can be 
used in assessing their performance as well as in improving them. The helium isoelectronic 
sequence provides such a class of systems, as also is the case for Hooke's atom sequence. 
However, as briefly pointed out above and as we shall demonstrate more carefully below, 
the former sequence is less general than what initially may be assumed, making it less in- 
formative for studies of approximate energy functionals. On the other hand, as we shall 



also demonstrate, Hooke's atoms provide a class of systems for which the members are "less 
similar" than is the case for the helium sequence, making a detailed study of this sequence 
highly relevant. 

We finally add that correlation energy is defined differently in wave function methods 
and in density functional theory. Since the numerical differences are not large and not will 
affect any of the conclusions of the present work, we shall not discuss this distinction further. 

2 The spherical two-electron systems 

The ground state of a spherically symmetric two particle system with the Hamiltonian 

^ = -\ (V? + V^) + v{n) + v{r,) + — , (3) 

I ^ ^ ri2 

in which the (local) single-particle potential v{r) is sufficiently attractive, is a bound singlet; 
hence, it satisfies pa = Pi3, where pa and p/j are the spin up and spin down densities, 
respectively; i.e., the spin polarization parameter ( = ^" /''^ vanishes. The spatial part of the 
corresponding Hartree-Fock wavefunction is of the form ^HF{f\,^) = 0(^1)0(^2), and the 
electron density is given in the form p{r) = 20(r)^ . It will be useful to note that |Vp(r)| = 
t {tf = i {if^ and VV = S + It Moreover, in the present context {^f = 
-^ {-£:) ■ For each of the density functionals considered below we provide the (sometimes 
considerably simplified) expression appropriate for the family of spherically symmetric two- 
particle systems presently studied. One reason for doing so is that the expressions for the 
various density functionals are scattered in the literature, in many instances the expression 
for a particular functional requiring consultation of several different sources whose notation 
is occasionally not fully consistent. The record is held by TPSS, the specification of which 
requires consultation of ref . [12], [33], [H], and [TT] . 
In our case we have 

{— for the helium sequence 
(4) 
^kr"^ for Hooke's atom. 



The correlation energy is given in terms of an expression of the form 

/■oo 

Jo 



where ec[p\ is a correlation energy density whose value at r may depend on the values of p, 

-£: as well as -^ at r. The dependence of the various functionals on the density is usually 

expressed via the Wigner-Seitz radius (the radius of a sphere containing one particle with 

uniform density p) 

( 3 \^/' 



\47rpy 
We follow the classification of ec[p] into 

• Zero-order approximation (OOA), depending on p only, 

• First-order approximation (lOA), depending on ^ as well, 

• Second-order approximation (20 A), depending also on -^. 

We note that this classification is closely related to the common classification into local 
density approximations (LDAs), generalized-gradient approximations (GGAs) and meta- 
GGAs. However, by restricting the classification according to the behavior for a spherical 
two-electron system, for a given functional we sometimes arrive at an expression that belongs 
to a lower class, in the classification specified above, than the general case. Our classification 
here applies to this reduced expression. 

Both as an additional check of the computations and to gain further insight into the 
behavior of the various density functionals considered, we examine the weak correlation limit 
(Z ^ oo or A; — > oo, for the Helium and Hooke's atom isoelectronic sequences, respectively) 
of the corresponding correlation energies. This can be done by using the asymptotic densities 

p(He) = exp(— 2Zr), Z ^ oo 

TT 



and 



/^\ 3/2 

p(Hooke) = 2 ( — j exp(— a;r^), lo = vk, A; ^ oo 



TT, 



In fact, the density, correct to first order in the interelectronic repulsion, can be written for 
the Hehum isoelectronic sequence by replacing the nuclear charge Z by the screened nuclear 
charge Zeg = Z — 5/16. Similarly, for Hooke's atom the variational principle yields for the 
trial function (^j exp(— ^r^) the equation uJIi^^\^J'^loJ^ —uj'^ = 0, where u = y/k. The 
solution, to first order, can be written in the form Ues ~ [V^— ny f ) • ^^ ia.ct, to study 
the weak correlation limit it is convenient to use the scaled radial coordinate f = Z^sr or 

~ 1/2 

r = uj^^ r, for the Helium and Hooke's atom isoelectronic sequences, respectively. In terms of 
these coordinates p{r)A'Kr'^dr becomes 8exp(— 2f)f^(if or -^ exp(— r^)f^(if, respectively. In 
terms of these scaled coordinates the expressions ^P"^''^? r dr*^"^^^ ^"^^ dr^P~^^^^ ^^^^ appear 
in several of the GGA and meta-GGA density functionals, obtain Z- (or k-) independent 
forms in the weak correlation limit. The weak correlation limit for the various density 
functionals examined below should be compared with the exact values jSHl 

E,(He) ^ -0.046663 + 2:2^ 

and 

00Q36Q 
Ee(Hooke) ^ -0.049703 + -^^^ • 

3 Zero order approximations 

3.1 Wigner's functional 

Wigner's first exploration of the correlation energy density as a means for estimating correla- 
tion energies in many-electron systems dates back to 1934 j^. Four years later he proposed 
the functional [^H] 

In the weak-correlation limit (i.e., when either Z ^* oo ox k ^ oo) Wigner's correla- 
tion energy for a two-electron system becomes — ^y|^ = —0.113, which is roughly twice 
as large as the corresponding (rather similar) limiting values of the helium and Hooke's 
atom isoelectronic sequences. A closer examination yields ii^^(He) ~ —0.113 + ^^- and 
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i?^(Hooke) ~ —0.113 + ^^^ji- These expressions are easily obtained analytically, noting that 
in this limit e^ ~ ~^^ + T^t^s- In the strong correlation limit ejf ~ —2^. 



3.2 The Gunnarsson-Lundqvist functional [39j 



e^^ = -0.0333 



l + x^)log(l + -) + :^-xJ-- 



(6) 



where Xp = j^. Since for r^ — ^ one has ef-^ ~ 0.0333 log(rs), it follows, using the 
asymptotic densities specified above, that 

^f^(He, Z ^ oo) ^ -2 ■ 0.0333 log(Z) - 0.09506 , 

and 

0333 

Ef ^(Hooke, k^oo)^ _uA|oo ^^^^^^ _ ^ ^^^ ^ 

The leading asymptotic terms for the helium and Hooke's sequences are related to one 
another by the substitution of k^ for Z, which is consistent with McWeeny's correspondence 
. In the low density (strong correlation) limit, r^ — ;> oo, hence e^^ ~ _a285_ 



3.3 The Brual-Rothstein functional 

This correlation energy density 

BR ^ 1 ^ 0-0289 

9.81 + 21.437p-i/3 0.284 + r, ^ ^ 

has the same form as Wigner's functional, eq. El The parameterization has been chosen to 

yield an exact correlation energy for the helium atom. 

Just like Wigner's functional, the weak correlation limit yields a common correlation 

energy for the hydrogenic and Hooke's sequences, E^^ = —0.204. 



3.4 The Perdew-Wang parameterization [II Ij 

This is the up-to-date local density approximation, that we write in the form 
er^W = -2co(l + «ir,)log 



1 + . 
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(8) 



where 

Co = 0.031091 

ai = 0.21370 

pi = 2coPi = 0.47232 

/32 = 2cop2 = 0.22308 

Ps = 2co/53 = 0.10187 

p4 = 2copi = 0.030652 . 

Pi, i = 1,2,3,4, (along with cq and ai) are given in ^T]. In the weak correlation limit 
{ts -^ 0) one obtains e™*-^(rs) ^ 2colog(/3iry^), hence 

Ef^(He) ^ -2colog(Z) - 0.05143 , 

and 

Ef ^(Hooke) ^ -J \og{k) - 0.07066 . 

These expressions agree with the computed correlation energies, presented in tables 1 and 
4, respectively. In the strong correlation limit, r^ — * oo, we obtain ejf"-*/ ^ _a434^ which is 
rather close to the corresponding limit of Wigner's correlation-energy density functional, eq. 



4 First-order approximations 

4.1 The Perdew86 functional [TUj 

Define 



and 



A = 0.002568 + 0.023266r, + 7.389 ■ lO^V,^ 



B = 1 + 8.723r, + 0.472rf + 0.07389rf . 
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Then C = 0.001667 + | and 

-0.0008129 dp „7/R 
C dr^ 

The correlation energy density for the homogeneous electron gas is parameterized as 
follows [Zj: For r^ < 1: 

e^°) = 0.0311 log(r,) - 0.048 + 0.0020r, log(r,) - 0.0116r, , 

[EJP\Re, Z ^oo) = -2- 0.0311 log(Z); ^^(Hooke, k ^ oo) = -^fii log(A;)] 



and for r, > 1: 



s 



m^ -0-1423 



l + 1.0529ryV0.3334r/ 
[Forr,^oo,eW^-^]. 

In terms of all the above, the correlation energy density is 

2 



6f^ = 4°)+Cexp(-$)(^y 



p-7/3 



In the Z ^ CO (or k -^ oo) limit the non local term vanishes, because in this limit $ ~ Z^^'^ 
($ ~ k^^^, respectively) so exp(— $) -^ 0. 

4.2 Becke's functional [8] 

Define B = (^)%"^/^. Then 

E, = -pi/3 (^ f^y + 0.00375 ■ 2^35(1 + 0.007 • 2'/'By'^'] . 
Now define z = — ^^. The expression for the correlation energy density is 



Ex 

er = -0.2pz^ (l - Ml±i) ^ . 

In the weak correlation limit it is found that 2; — > 0, hence ef ~ — O.lp^;^. Numeri- 
cal integration over the asymptotic densities yields El^^^{B.e, Z —>■ 00) = —0.07445 and 
Ef ^^(Hooke, k^oo) = -0.22176. 
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4.3 The Wilson-Levy functional [32] 
This functional, 



ML 



-0.74860 + 0.06001 



P 



-4/3 



p-4/3 + 



(9) 



3.60073 + 2.2679 

is a non-local generalization of Wigner's functional. Numerical integration over the asymp- 
totic densities yields Ef ^(He, Z ^ oo) = -0.04803 and E^^^(Hooke, k ^ oo) = -0.07817. 



4.4 The Perdew, Burke and Ernzerhof functional [14j 

For a spherically symmetric singlet state this correlation-energy density is of the form 



e— (r.) + i/(r„t) 



(10) 



e""'*-'"(rs) is given by eq. |Hl and 



H{rs,t)=-f log {l + 6f 



1 + Af 



1/6 



where t = i (|) ^ P^^^^f and A 



exp 



l + At^+ (At2)2 
1-1 



unif 



-1 



Here, 7 = Cq (c/. eq. |HI) 



and 5 = 2.1461. The Z ^ oo logarithmic leading term of e""*-^(rs) is canceled by an 
equal term with opposite sign in H{rs,t). A recent comparative study ^ suggests that this 
correlation energy density functional is superior to all other density functionals considered 
in that study. The weak correlation limit requires numerical integration over the asymptotic 
densities, yielding E^^^{B.e, Z -^ cxd) = —0.04789, in agreement with [T3j and [TTj, and 
Ef ^^(Hooke, k-^oo) = -0.08131. 



4.5 The Filatov-Thiel density functional 

Filatov and Thiel proposed two closely related functionals, that we shall designate FT97 (cf. 
53^]) and FT98 (cf. |13]). Whereas FT97 belongs in the present subsection, FT98 contains 
a dependence on the Laplacian of the density and belongs in the following subsection. 
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FT97 "[...] is based on a model of the Coulomb hole distribution function and accurately 
describes the correlation energy of the uniform electron gas and the atoms He to Ar" (44j . 
Define 



and let 



Let 



and 



'97 



1 + 1.622199yo + (0.489958yo) 
1 + 1.094530i/o/r, 

0.02940rs 



2 



■exp[-(0.489958i/o) ] • (H) 



/^ 



{[l + 0.370365 exp(-0.10018rf^)] Fgy}^ 



/(/i) = -exp(/i)Ei(/x) 

where -E'i(x) is the exponential integral. 

The FT97 correlation density functional |^ El] is 



e™ = 0.015545 |/(/.) + ^^^^!^[/^/(^) + 1]} • (12) 

In the weak correlation limit E^^^^CRe) ^ -0.05648 + 2^, and Ef"^(Hooke) ^ 

0.07 



-0.07076 + ^°^ 



4.6 The One-parameter progressive CS-type functional [46| 

This correlation energy density is a variant of the Cole-Salvetti formula, discussed in sub- 
section 5.1. Defining x„ = 2^^^p~^^^-£:, this functional can be written in the form 

OP noonoK /3„/3 + 0.3789 , , 

'= =-''-™^«''/3^,+ 1.1284/3J, + 0.3183«, <"' 

where 

3\^/3 0.0084x2 



/?„^ = 0.93934p 



1/3 



3 — 



Air J l + 0.0252x^log(a;^ + y'x2 +1)_ 

In the weak correlation limit ^^^^(He) ^ -0.08511+^, and£;f^(Hooke) ^ -0.09504+ 



0.074 



14 



4.7 Perdew, Kurth, Zupan and Blaha functional [35] 

This functional depends on the kinetic energy density, and is therefore a meta-GGA func- 
tional. However, for the ground state of the two-electron system it reduces to 

ef^-^ = 1.53(6^^ - eT) , (14) 

where e^^^ is given by eq. ^1 and 

Here g^"-*/'"P has the same form as ejf"*-^, eq. |Hl with the parameters Cq, ai and /3j, i = 1,2,3,4 
being replaced by 

Co = 0.015545 

ai = ai2^/^ = 0.25889 

;5i = 2co/?i"^2^/^ = 0.49273 

^2 = 2co/372^/3 = 0.24278 

'^s = 2co/?3^2i/2 _ 0.14801 

^4 = 2co/?r2^/^ = 0.03085 , 

respectively, cq, ai and /Jf^, i = 1,2,3,4 are presented in ref. [TT], Table I. The fractional 
powers of 2 take care of the fact that the spin up density is a half of the total density, so r^ 
should otherwise have been replaced by rs2^^^. 

H'^^{rs,t) is obtained from H{rs,t) by multiplying it by 1, replacing A by 



A''P = S 



exp — 1 

V 7 / 



H -1 



and replacing t by t"^ = t\^. 

Inspecting the reparametrized expression we note that 



^unif,up/ \ ^ _^unif I 
15 



This is not apparent in terms of the parameterization given in JT] . 

Numerical integration over the asymptotic density functionals yields E^^^^{B.e, Z — > 
oo) = —0.05525, in agreement with the value reported in the erratum to ref. |17,, and 
Ef ^^^(Hooke, k^ oo) = -0.08468. 



4.8 The TPSS functional 

This functional is of the form 

where e^^^^ is given by eq. El The asymptotic correlation energies are Ej^^^{B.e, Z — »• 
cxd) = —0.04901, in agreement with the value reported in the erratum to ref. ^7], and 
Ej^^^(Hooke, k^oo) = -0.06843. 

5 Second-order approximations 

5.1 Colle-Salvetti formula 

The Colle-Salvetti formula [5_^ for the correlation energy depends on the Hartree-Fock two- 
particle density matrix. For the ground state of a two-electron atom it can be written in 
terms of the Hartree-Fock density, providing a correlation energy density 

.f^ = -0.04918 ' + 0-132p-V^exp(-0,2533p-'/») r^ 

l + 0.349p-i/3 ' ^ ^ 



where 

1 



£p 2 dp 1 fdp^ 
dr^ r dr p \dr ^ 



Integration over the asymptotic densities, which can be carried out analytically, yields 
Ef ^(He, Z ^oo) = -0.01941 and ^^(Hooke, k -^ oo) = -0.001772, respectively. 
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5.2 Lee, Yang and Parr functional [3] 



This correlation-energy density is obtained by replacing tcs, in eq. ^1 by 

p/, Y / d^p 'ldp\ iY / dp \ 



^ f(Pp 2dp' 
24 \ dr"^ r dr 



ji(3^2)2/3^5/3_^^,_^^__^,_ 



72p \dr ) 

For the Z ^ oo (hydrogenic) density the LYP correlation energy functional can be 
integrated analytically, to yield the correlation energy -0.0566888. This value, obtained in 
ref. [TTj, appears to be too high (in absolute value) to be consistent with the values presented 
in Table 3. Closer examination reveals that the asymptotic behavior of the LYP correlation 
energy for the helium isoelectronic sequence is given by i?c(He) ~ —0.0566888 + 0.009 °^^^ . 
This expression, which is consistent with the numerical results, can be obtained analytically 

. It is perhaps worth noting that using the screened hydrogenic density, in which Z is 



replaced by Z — ^, the LYP functional yields, for Z = 2, the correlation energy -0.04388, 
compared to the value -0.04378 obtained with the Hartree-Fock density. For Hooke's atom 
isoelectronic sequence the asymptotic behavior of the correlation energy, at A; — >■ cxd, is found 



to be Ec(Hooke) ^ -0.035501 + 0.00012 



A:i/4 • 



5.3 FT98 [45] 

This is a modified version of the FT97 functional, extended to include the Laplacian of the 
density. It can be written by replacing F97, eq. [TTl by 

l + (0.083726y„^2 



1 + 0.05158l2/^/r, 



exp[-(0.083726i/^) 



where 



Vcj 



iP 



-8/3 



' dp^ 



+ 



0.89672 



^187r)2/3^ \dr) (187r)4/3 

In the weak correlation limit ii^J^^^(He) 
0.22 

fcl/4- 



P 



-10/3 



1 (dp" 

P \dr ^ 



'd^p 2 dp" 
dr^ r dr 



-0.08173 + ^, and £:™(Hooke) 



-0.08706 + ^22 
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6 Performance of the various density functionals 

The Hartree-Fock densities were evaluated for both the hehum and Hooke's isoelectronic 
sequences, using the procedure described in ref. jlH]- These densities were used to compute 
the correlation energies in terms of the various density functionals specified above. The 
results are presented in Tables 1-6, where comparisons with previous pertinent results are 
documented, as well as in Figures 1-8. When comparing the different results it is important 
to pay attention to the ranges of the correlation energy scales in the different Figures. Thus, 
Fig. 1, presenting the OOA correlation energies for the helium isoelectronic sequence, has a 
range five times as large as Fig. 2 (crude lOA), whose range is three times as large as that 
of Fig. 3 (better lOA). The 20A results, presented in Fig. 4, do not improve on the best 
lOA results in Fig. 3. The overall trend is similar for Hooke's atom isoelectronic sequence, 
but much less pronounced. Thus, the range of Fig. 5 (OOA) is twice that of Fig. 6 (crude 
lOA), where the range of deviations of the dft correlation energies from the exact values is 
only a factor of 1.5 broader than that in Fig. 7 (better lOA). Again, the 20A correlation 
energies do not improve on the best lOA results. 

The four OOA functionals presented in Fig. 1 for the helium isoelectronic sequence fail 
to provide an adequate estimate of the correlation energy. Only BR provides the correct 
correlation energy for the helium atom, because it had been calibrated to do so. Wigner's 
is the only OOA- type functional providing a mild dependence on the nuclear charge (Z), 
that mimics the exact correlation energy, though it is more than twice as large. In fact, 
the Wigner and BR correlation-energy functionals are the only local density functionals 
considered, that yield finite correlation energies for Z ^ oo. Needless to say, since PW is 
the definitive LDA limit, its failings should be taken as a genuine assessment of this level of 
approximation. That other OOA-type correlation energies are in closer agreement with the 
exact results should be interpreted as fortuitous. It is regrettable that a "first-principles" 
benchmark lOA correlation energy density cannot be specified in a similar way, or has not 
so far been proposed. For Hooke's atom all zero order approximations considered, presented 
in Fig. 5, are roughly equally poor. At intermediate values of the force constant the ordering 

18 



of the various approximations relative to the exact correlation energy is similar to that for 
the helium isoelectronic sequence. 

The first-order functionals were split into two classes according to their behavior for the 
helium isoelectronic sequence. The more poorly behaved functionals, P86, OP and Becke, 
are presented in Fig. 2. They exhibit a strong dependence on the nuclear charge. The five 
functionals presented in Fig. 3 all have a reasonable slope, only mildly steeper than that 
of the exact correlation energy. Both PBE and Becke have been parameterized to yield 
the exact correlation energy for the helium atom. For Hooke's atom isoelectronic sequence 
most lOA results, except P86, cluster together, being roughly twice as large as the exact 
correlation energy. We note that TPSS and FT97 are remarkably close to one another; so are 
PBE and PKZB. Neither pair remains close for the helium isoelectronic sequence. Rather, 
PBE and TPSS behave similarly; so do FT97 and PKZB, even more closely but further away 
from the exact correlation energy. These contradictory patterns serve as a warning against 
assigning excessive significance to limited scope examination of the comparative performance 
of different density functionals. 

The three second-order approximations considered are presented, for the helium isoelec- 
tronic sequence, in Fig. 4. LYP has the overall most satisfactory behavior, although for 
Z < 6 CS is in closer agreement with the exact correlation energies. For Hooke's sequence 
(c/. Fig. 8), none of the second-order approximations provides a good agreement with the 
exact correlation energies. LYP and CS, which are conceptually related, are remarkably close 
to one another, along Hooke's atom isoelectronic sequence; however, this is not true for the 
helium isoelectronic sequence. 

As already pointed out, Hooke's atom isoelectronic sequence allows the examination of 
the weak binding (high correlation) limit. A; — ^ 0. The computed correlation energies for the 
range of lower force constants suggest that 



JT'exact 



l.lk^ 
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and 

^^- ^ l + 2.18ks , 

Tpexact ' 

although, in both cases, the correlation energies are considerably higher than the exact 
correlation energies, even for the lowest force constants considered. These results may be 
taken to suggest that at very low force constants the local density approximation becomes 
exact, since the density is spread out and its gradient becomes negligibly low. A similar 
behavior is observed for the ratios pl^aa , pexad and -nfes^. This feature has not been 
deliberately built into any of these approximations. It is, therefore, noteworthy that it holds 
to the extent that it does. 

Our study clearly suggests that the approximate density functionals in general perform 
better for the helium sequence than for Hooke's atoms. In order to understand this behavior 
we show in Fig. 9 scaled electron densities, i.e., p/c"^ as well as ^/c^ and ^/c^, as functions 
of X = c ■ r with c being Z — 5/16 for the He sequence and k^^^ — 0.06649 for Hooke's atoms. 
These scaling parameters were derived at the bottom of section 2, above. 

A number of observations can be made from the figure. First, the curves for the He 
sequence are all remarkably similar, except for the region r ^ 0, where the cusp problem 
becomes visible, in particular for the second-order derivative. An analogous similarity is not 
found for Hooke's atoms. The curves corresponding to the highest values of k are indeed very 
close to the expected k ^ oo limit, but the range of k is broad enough to allow significant 
deviations from this limiting behavior. It follows that the He sequence does not provide 
a useful testing ground for exploring approximate density functionals, as long as it is not 
extended with further test systems. 

Second, the electron density is more confined in space for the He sequence than for 
Hooke's atoms. For the former, the cusp problem may be important since the region for 
r -^ is one of a large energy density although its volume is vanishing. A similar region 
does not exist for Hooke's atoms, i.e., for those systems the electron density is not maximized 
in a region of diverging potential. This suggests a further reason for using Hooke's atoms 
as model systems for studying properties of electron densities in regions where the potential 
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stays finite. This is exactly what is relevant when developing approximate density functionals 
aimed at describing the properties of chemical bonds. 

7 Conclusions 

A word is appropriate about the choice of density functionals investigated. We did try to 
be fairly exhaustive. Some functionals were included primarily to respect the pioneers, most 
notably Wigner, although we have not made justice to contributions by Pines jSU], Gell- 
Mann and Brueckner jM], Hedin [521, ^^^ many others. We did exclude a couple of recent 
functionals either because the correlation energies they provide for helium were judged to 
be poor, because no result was available to verify our implementation, or because they were 
considerably harder to implement. In the cases excluded we did not have a compelling 
reason to expect a particularly good performance for the helium sequence, let alone Hooke's 
sequence. We may have excluded highly worth-while functionals by oversight or ignorance, 
for which we apologize. 

It has already been noted by several authors that most density functionals fail to yield an 
adequate dependence of the correlation energy on the nuclear charge, Z, along isoelectronic 
sequences. However, some, such as PBE, TPSS, CS and LYP are not too bad. None of 
that is true for Hooke's sequence, suggesting that the ultimate aim of developing a universal 
density functional has not yet been attained, not even for spherically-symmetric two-fermion 
systems. By analysing the electron densities for the two classes of systems, we found that 
the He sequence does not provide a sufficiently general class of systems for assessing or 
improving the quality of density functionals, whereas this is to a larger extent the case for 
Hooke's atom isoelectronic sequence. A feasible approach to improving the universality of 
correlation energy density functionals could be to examine their reparametrization, allowing 
for a dependence on the one-body potential. Since this potential is uniquely determined by 
the ground-state Hartree-Fock density, this approach is still consistent with the Hohenberg- 
Kohn framework of dft. 
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z 


Tpexact 


EPW 


E^"- 


EBR 


EW 


2 


-0.042044 


-0.1 125 [a] 


-0.1385 


-0.04204 


-0.0959 


3 


-0.043498 


-0.1346[b] 


-0.1666 


-0.05922 


-0.1015 


4 


-0.044267 


-0.1 504 [a] 


-0.1862 


-0.07276 


-0.1043 


5 


-0.044736 


-0.1628 


-0.2037 


-0.08376 


-0.1059 


6 


-0.045054 


-0.1730 


-0.2137 


-0.09290 


-0.1071 


7 


-0.045281 


-0.1817 


-0.2240 


-0.10065 


-0.1079 


8 


-0.045452 


-0.1893 


-0.2330 


-0.10730 


-0.1085 


9 


-0.045586 


-0.1960 


-0.2409 


-0.11309 


-0.1090 


10 


-0.045692 


-0.2020[a] 


-0.2499 


-0.11819 


-0.1093 


20 




-0.2424 


-0.2944 


-0.14834 


-0.1110 


50 




-0.2971 


-0.3555 


-0.17618 


-0.1121 



[a] ref. [I2|. [b] ref. 



Table 1: Zero-order-approximation correlation energies for the Helium isoelectronic sequence. 
The approximate density functionals are defined in the text. 
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to 

00 



z 


rpexact 


EPBE 


pPKZB 
^c 


EOP 


pFT97 
^c 


^TPSS 


pB88 
^c 


^P86 


jp^WL 


2 


-0.042044 


-0.04202 [b] 


-0.04727[b] 


-0.05178[c] 


-0.04623 [d] 


-0.04306 [b] 


-0.04189[e] 


-0.04390 [f] 


-0.04204[a] 


3 


-0.043498 


-0.04478 [b] 


-0.05063 [b] 


-0.06089 


-0.04952 


-0.04574[b] 


-0.04985 [e] 


-0.04542 [f] 


-0.04418 [e] 


4 


-0.044267 


-0.04606 [b] 


-0.05224[b] 


-0.06607 


-0.05119 


-0.04698 [b] 


-0.05462 [e] 


-0.04868 [f] 


-0.04520 [e] 


5 


-0.044736 


-0.04676 


-0.05316 


-0.06943 


-0.05222 


-0.04767 


-0.05782 


-0.05324 


-0.04579 


6 


-0.045054 


-0.04720 


-0.05374 


-0.07178 


-0.05288 


-0.04811 


-0.06012 


-0.05865 


-0.04618 


7 


-0.045281 


-0.04748 


-0.05441 


-0.07351 


-0.05338 


-0.04839 


-0.06186 


-0.06457 


-0.04645 


8 


-0.045452 


-0.04768 


-0.05441 


-0.07485 


-0.05375 


-0.04860 


-0.06322 


-0.07079 


-0.04665 


9 


-0.045586 


-0.04782 


-0.05463 


-0.07591 


-0.05404 


-0.04874 


-0.06431 


-0.07717 


-0.04681 


10 


-0.045692 


-0.04793 


-0.05478 


-0.07677 


-0.05428 


-0.04885 


-0.06521 


-0.08360 


-0.04693 [a] 


20 




-0.04826 


-0.05536 


-0.08080 


-0.05535 


-0.04923 


-0.0695 


-0.14274 


-0.04749 


50 




-0.04824 


-0.05550 


-0.08335 


-0.05602 


-0.04928 


-0.07240 


-0.24874 


-0.04781 



[a] ref. iJjj. [b] ref. jHSj. [c] ref. [46j. [d] ref. |1H]. [e] ref. ^. [f] ref. |lil|. 



Table 2: First-order-approximation correlation energies for the Helium isoelectronic sequence. The approximate density 
functionals are defined in the text. 



z 


rpexact 


Ef^ 


ELYP 


^FT98 


2 


-0.042044 


-0.04158[c] 


-0.04378 [a] 


-0. 04643 [d] 


3 


-0.043498 


-0.0439 [g] 


-0. 04755 [e] 


-0.05254 


4 


-0.044267 


-0.04421 [g] 


-0.04905[a] 


-0.05804 


5 


-0.044736 


-0.04385 


-0.04972 


-0.06122 


6 


-0.045054 


-0.04326 


-0.05004 


-0.06360 


7 


-0.045281 


-0.04260 


-0.05018 


-0.06546 


8 


-0.045452 


-0.04191 


-0.05022 


-0.06697 


9 


-0.045586 


-0.04124 


-0.05023 


-0.06821 


10 


-0.045692 


-0.04060 


-0.05021 [a] 


-0.06925 


20 




-0.03591 


-0.04992 


-0.07464 


50 




-0.03001 


-0.05008 


-0.07862 



[a] ref. [121 • [c] ref. HSj. [d] ref. ^ [e] ref. 



ref. m. 



Table 3: Second-order-approximation correlation energies for the Helium isoelectronic se- 
quence. The approximate density functionals are defined in the text. 
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k 


rpexact 


EPW 


E^^ 


EBR 


E^ 


0.000300891 


-0.01691 


-0.02865 


-0.0303 


-0.00392 


-0.0388 


0.00133497 


-0.02198 


-0.03857 


-0.0425 


-0.00609 


-0.0506 


0.01 


-0.02904 


-0.05478 [h] 


-0.0633 


-0.01071 


-0.0665 


0.25 


-0.03843 


-0.08614[i] 


-0.1045 


-0.02424 


-0.0871 


1.0 


-0.04134 


-0.1013 [h] 


-0.1243 


-0.03321 


-0.0935 


16. 


-0.04527 


-0.1342 


-0.1663 


-0.05784 


-0.1023 


100. 


-0.04684 


-0.1575 [h] 


-0.1952 


-0.07847 


-0.1060 


10000. 


-0.04877 


-0.2208 [h] 


-0.2699 


-0.13417 


-0.1105 



[h] ref. t32j. [i] ref. ^ 



Table 4: Zero-order-approxiination correlation energies for Hooke's atom isoelectronic se- 
quence. The approximate density functionals are defined in the text. 
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w 



k 


rpexact 


j^PBE 


j^PKZB 


E^P 


j^FT97 


j^TPSS 


jj^BSS 


^P86 


j^WL 


0.000300891 


-0.01691 


-0.02148 


-0.01895 


-0.00970 


-0.01904 


-0.01835 


-0.00860 


-0.02508 


-0.0385 


0.00133497 


-0.02198 


-0.02782 


-0.02525 


-0.01446 


-0.02483 


-0.02416 


-0.01247 


-0.03244 


-0.0475 


0.01 


-0.02904 


-0.03705 


-0.03470 


-0.02341 


-0.03292 


-0.03260 


-0.01951 


-0.04216 


-0.0577 


0.25 


-0.03843 


-0.05118[j] 


-0.04968 


-0.04265 


-0.04484 


-0.04523 


-0.03448 


-0.05434 


-0.0684 


1.0 


-0.04134 


-0.05748 


-0.05546 


-0.05159 


-0.04924 


-0.04981 


-0.04168 


-0.05833 


-0.0712 


16. 


-0.04527 


-0.06520 


-0.06522 


-0.06768 


-0.05652 


-0.05708 


-0.05547 


-0.06712 


-0.0746 


100. 


-0.04684 


-0.06960 


-0.07028 


-0.07594 


-0.06023 


-0.06057 


-0.06315 


-0.07683 


-0.0759 


10000. 


-0.04877 


-0.07659 [k] 


-0.07861 


-0.08816 


-0.06624 


-0.06569 


-0.07563 


-0.12609 


-0.0774 



[j] ref. |2H]. [k] ref. ^. 



Table 5: First-order-approximation correlation energies for Hooke's atom isoelectronic sequence. The approximate 
density functionals are defined in the text. 



k 


Tpexact 


E^' 


ELYP 


E™ 


0.000300891 


-0.01691 


-0.01068 


-0.01069 


-0.01141 


0.00133497 


-0.02198 


-0.01559 


-0.01561 


-0.01563 


0.01 


-0.02904 


-0.02367 


-0.02369 


-0.02239 


0.25 


-0.03843 


-0.03542 


-0.03502[m] 


-0.03488 


1.0 


-0.04134 


-0.03837 


-0.03762 


-0.04057 


16. 


-0.04527 


-0.03913 


-0.03799 


-0.05183 


100. 


-0.04684 


-0.03636 


-0.03573 


-0.05879 


10000. 


-0.04877 


-0.02397 


-0.02871 [m] 


-0.07296 



[in] ref. 



Table 6: Second-order-approxiination correlation energies for Hooke's atom isoelectronic 
sequence. The approximate density functionals are defined in the text. 
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Figure 1: Zero-order-approximation correlation energies for the helium isoelectronic se- 
quence. The approximate density functionals are defined in the text. 
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Figure 2: First-order- approximation correlation energies for the helium isoelectronic sequence 
for less-well performing functionals. The approximate density functionals are defined in the 

text. 
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Figure 3: First-order- approximation correlation energies for the helium isoelectronic sequence 
for well performing functionals. The approximate density functionals are defined in the text. 
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Figure 4: Second-order-approximation correlation energies for the helium isoelectronic se- 
quence. The approximate density functionals are defined in the text. 
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Figure 5: Zero-order-approximation correlation energies for Hooke's isoelectronic sequence. 
The approximate density functionals are defined in the text. 
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Figure 6: First-order-approximation correlation energies for Hooke's isoelectronic sequence 
for less-well performing functionals. The approximate density functionals are defined in the 

text. 
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Figure 7: First-order-approximation correlation energies for Hooke's isoelectronic sequence 
for well performing functionals. The approximate density functionals are defined in the text. 
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Figure 8: Second-order-approximation correlation energies for Hooke's isoelectronic se- 
quence. The approximate density functionals are defined in the text. 
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Figure 9: The electron density and its derivatives for the hehum isoelectronic sequence (left 
panels) and Hooke's atom isoelectronic sequence (right panels). All curves are given as 
functions of a scaled radial coordinate x = c- r, with c being Z — 5/16 for the He sequence 
and k^^^ — 0.06649 for Hooke's atoms. For Hooke's atoms, the uppermost curves for small x 
in the upper panel and the lowest curves for small x in the other two panels correspond to 
the largest values of k. 
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